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PHENOMENOLOGTCAL  MODELLING  OF  P0LYWELL‘“/SCIF 
MULTI-CUSP  INERTIAIr-ELECTROSTATIC  CONFINEMENT  SYSTEMS 

L  INTRODUCTION  AND  BACKGROUND 


In  an  earlier  paper^  (hereinafter  referred  to  as  ”[1]"))  &  comprehensive  study  was 
made  of  the  recirculation  and  losses  of  electrons  in  their  flow  through  simple  inverse 
power— law  potential  wells  bounded  by  similarly  inverse  power-law  dependent  magnetic 
fields.  This  study  examined  electron  flow  and  loss  behavior  in  the  simplest  approximation 
invoked  to  describe  Polywell*“  confinement  systems.  The  importance  of  this  study,  and  of 
thd  present  paper,  is  that  the  power  balance  in  Polywell^’”  systems  is  determined  entirely 
by  the  rate  of  electron  losses;  if  these  are  large,  then  the  system  can  not  yield  net  power. 

As  noted  in  [I],  a  large  body  of  work  has  been  undertaken  over  the  past  35  or  so 
years  in  the  study  of  general  cusp  confinement  of  plasmas^,  s.  s.  Nearly  all  of  this  has 
examined  single  particle  electron  (or  ion)  motion  or  the  motion  of  particles  in  neutral 
plasmas  within  cusped  magnetic  systems,  generally  without  internal  electric  potential 

i 

fidds.  Furthermore,  almost  all  of  this  work  on  plasmas  has  been  limited  to  plasmas  in 
local  thermodynamic  equilibrium  (LTE);  none  applied  to  non-neutral  systems.  Nearly  all 
were  without  internal  E  fields;  all  particles  studied  were  at  constant  (fixed)  total  kinetic 
energy  throughout  the  system.  And,  almost  all  of  the  prior  work  focussed  on 
spindle-shaped  biconic  cu8p8>  — flbich  are  uniqudy  unable  to  satisfy  the  configuration 
criteria  necessary  for  effective  electron  confinement  in  Polywell*”*— type  multicusp  systems. 

These  spedal  polyhedral  configurations^^  *  allow  formation  of  stable  deep 
dectrostatic  potential  wdls  in  dectron/ion  mixtures  by  radial  injection  of  energetic 
dectrons  into  plasmas  with  net  negative  charge.  Thdr  spatial  characteristics  combined 
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with  the  distribution  of  the  internal  E  field  arising  from  electron  injection,  and  the  induced 
predominantly— radial  ion  motion,  can  result  in  long  confinement  times  (due  to  many 
recirculation  passes  through  the  well  before  escape)  for  the  dectrons^.  This  is  the  sine  qua 
non  for  minimum  electron  energy  losses,  under  the  conditions  of  large  ion  energy  and  core 
density  achieved  by  the  trapping  of  ion  motion  by  the  electron— driven  internal  potential. 

In  this  note  the  problem  of  electron  motion,  confinement  and  losses  is  analyzed  in 
the  non— LTE,  electron— rich  plasma  system  with  anisotropic  radial  energy  in  both  species, 
in  multicusp  Polywell*™  geometry.  The  bulk  problem  is  treated  as  one— dimensional,  with 
arbitrary  spatial  indices  of  radial  B  and  E  field  variation.  Here,  unlike  [I],  the  B  field  is 
taken  to  have  a  "rollover"  shape,  consistent  with  the  results  of  computer  calculations  of 
actual  field  distributions  in  polyhedral  systems,  that  follows  a  simple  power— law  within  the 
system  but  departs  from  it  as  the  system  boundary  is  approached  from  within. 

An  important  feature  of  the  Poljrwell*™  system  is  the  possibility  of  inducing  B  field 
distortion  by  bulk  diamagnetic  collective  effects  arising  from  diamagnetic  currents  induced 
in  the  motion  of  electrons  by  their  interaction  with  the  basic  field  geometry.  These 
diamagnetic  currents  will  act  to  exclude  the  B  field  from  the  interior  of  the  system;  an 
effect  which  will  lead  to  a  B  field  shell,  compressed  to  a  ^  =  1  condition,  surrounding  and 
confining  the  central  electron  cloud  in  the  system.  These  effects  are  crudely  modelled  as 
they  influence  this  1-D  spatial  variation  of  B  field,  but  electron  entry  into  single  cusp 
volumes  includes  elements  of  tit  geometric  effects  of  the  real  3-D  configuration.  Electron 
motions  in  mirror  reflection  regions  are  analyzed  on  the  usual  basis  of  conservation  of 
magnetic  moment  of  the  electron  at  entry  into  the  confining  cusp.  However,  turning  points 
of  this  motion  are  modified  to  account  for  the  effects  of  diamagnetic  currents  on  this 
process. 
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In  the  discussion  here,  it  is  helpful  to  define  several  terms  for  critical  radii  within 
the  system.  The  electrons  are  moddled  as  ori^nating  at  the  system  center,  at  r  =  0, 
moving  through  a  non-iadiabatic  region  to  that  radius  (r^^)  at  which  the  electron  "/9  =1", 
thence  to  that  radius  (r^)  at  which  adiabatic  reflection  will  occur  in  the  cusp  fidds,  and 
finally  to  the  outer  radius,  r  =  R,  of  the  system. 

n.  MIRROR/WIFFLE-BALL  MODE  OPERATION 

There  are  two  distinct  modes  of  operation  in  the  system.  The  first  operating  scheme 
is  to  establish  the  system  with  maximum  6  field,  and  turn  on  a  small  electron  current  at 
maximum  (desired)  electron  energy,  followed  by  injection  of  ions  controlled  so  as  to  yield 
the  desired  central  virtual  anode  height  in  the  system.  It  has  been  shown  by  Bussard  and 
King^o  that  this  control  must  be  relatively  precise  (e.g.  within  a  &ctor  of  2— 4x)  in  order  to 
avoid  "blowout"  of  the  well  by  excessive  ion  density.  Given  a  controllable  ion  source,  and 
similarly  controllable  electron  injection,  the  system  can  be  driven  to  higher  core  density  by 
injection  of  ever-increasing  electron  current  with  corresponding  ion  input. 

This  mode  of  operation  ensures  that  the  system  will  start  as  a  mirror  confinement 
system  for,  at  low  electron  density,  <r|,>  =  *8  very  small  compared  to 

<r^>  =  (r^/R)  and  the  only  physics  initially  able  to  retain  the  electrons  is  that  of 
reflection  in  the  cusp  fidds;  no  diamagnetic  0^1  fidds  have  yet  been  established.  Once 
started  in  this  mirror  refiectiofli^R)  mode,  continued  increase  in  dectron  injection 
current  drives  the  0=1  radius,  <r^>,  to  larger  values  so  that  some  internal  reflection 
within  the  diamagnetic  sphere  adds  to  the  confinement  due  to  continuing  MR  effects.  It  is 
important  to  note  that  the  dectron  redrculation  due  to  the  MR  contribution  is  severdy 
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limited  by  the  considerable  depth  of  the  negative  potential  well  created  by  the  (desired) 
excess  electron  charge  in  the  system. 

In  contrast  with  electron  motion  under  conservation  of  magnetic  moment  in 
fidd--free  cusp  systems,  the  confining  effect  of  the  magnetic  moment  conservation  process 
is  greatly  reduced  here  by  the  presence  of  the  internal  E  field.  This  acts  to  accelerate  the 
electrons  outwards  in  the  system,  and  adds  radial  energy  to  their  motion  almost  as  fast  as 
the  cusp— trapped  transverse  energy  increases  and  takes  it  from  radial  motion.  In  fact,  if 
the  well  depth  were  to  become  exactly  equal  to  the  electron  injection  energy  (a  physical 
impossibility  so  long  as  ions  are  present),  the  cusp  fields  would  not  confine  the  electrons  at 
all.  The  MR  mode  works,  but  only  poorly  if  the  system  has  a  deep  potential  well. 

The  general  model  of  electron  motion  in  the  system  for  this  mode  of  startup  and 
operation  was  discussed  and  illustrated  further  in  [I],  which  showed  the  linear  ap¬ 
proximation  used  for  variation  of  the  adiabatidty  radius,  <r^>,  in  modelling  the 
transition  from  the  MR  mode  to  the  diamagnetic  0=1  mode.  Examining  this  mode  in 
some  more  detail,  once  the  electrons  have  reached  a  radial  distance  (r^)  at  which  their  local 
gyro-magnetic  radius  is  sufGdently  small  that  the  total  area  of  all  of  their  gyro  radius 
drdes  from  each  of  the  cusps  is  equal  to  the  area  of  the  spherical  surface  at  that  radius, 
they  will  enter  one  or  another  cusp  in  their  outward  motion.  If  r,  is  less  than  this 
"transition  radius"  r^  the  electrons  move  non-adiabatically  without  confinement,  and  enter 
the  spatial  domain  of  each  cus^Rt  r^.  No  diamagnetic  effects  are  yet  functional  at  this 
point,  as  the  gyro  radius  of  electrons  within  r^  <  r^  is  greater  than  the  gyro  radius  at  r^, 
and  coherent  currents  counter  to  the  external  B  fields  can  not  be  generated  and  sustained. 

Once  r^^  becomes  greater  than  r^  the  electrons  inside  r^^  will  "see"  a  lossy  sphere 
with  surface  loss  area  (due  to  gyro  radius  area  "holes")  that  is  smaller  than  the  sphere  area 
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at  and  they  will  be  confined  by  internal  reflection  within  this  diamagnetic  ^  =  1  sphere. 
The  average  number  of  transits  an  electron  will  make  within  this  sphere  at  r^  is  just  the 
ratio  of  the  sphere  area  to  the  total  loss  area.  As  r^  grows  larger  (with  increasing  electron 
current  input,  and  increasing  central  density  —  with  requisite  additional  ions  to  maintain 
charge  balance)  this  confinement  grows  larger,  because  the  fractional  hole  loss  area 
becomes  smaller  due  to  increasing  B  field  strength  at  larger  radial  positions.  This  can  be 
thought  of  as  the  confinement  of  a  particle  inside  a  perfectly-reflecting  spherical  shell 
perforated  by  holes  —  like  the  "wiffle  ball”  toy;  and  it  is  called  "wiffle— ball”  (WB) 
confinement.  The  average  number  of  transits  within  r^^  is  defined  here  as 

At  some  radius  (r^)  set  by  the  mirror  adiabaticity  requirement  that  the  local  B 
field  change  only  little  over  a  local  gyro  radius,  the  electron  is  trapped  in 
’’mirror— reflection”  (MR)  oscillation  within  the  cusp.  Electron  trapping  in  this  (MR) 
mode  is  by  the  usual  form  of  mirror  reflection  coefficient  for  motion  in  a  single  cusp  mirror 
system^  .  This  can  be  used  to  find  a  measure  of  the  average  number  of  transits  G.  that 
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an  electron  will  make  between  r^  <  r  <  R  within  this  cusp  mirror  geometry. 


The  form  for  mirror  reflection  coefficient  involves  a  term  which  is  the  square  of  the  sine  of  the 

minimum  trapping  angle  (relative  to  the  B  field)  of  electron  motion  in  the  B  field  at  r^,  the  trapping 

radius,  for  their  escape  at  radius  r4PB.  In  the  case  of  a  system  with  electrons  at  constant  energy  and  a 

B  field  varying  as  <r>*”  this  term  varies  as  present  case,  the  simultaneous  radial 

variation  of  the  E  field  reduces  the  reflection  coefficient  by  virtue  of  the  acceleration  of  electrons  out  of 

the  system  by  the  potential  well  gradient.  If  this  E  field  is  equal  in  magnitude  to,  and  varies  with  the 

same  functional  form  as  the  B  field,  no  mirror  trapping  will  be  possible  in  the  system.  If  the  potential 

well  is  less  deep  than  the  electron  injection  energy,  trapping  is  still  powible,  but  with  the  reflection 

coefficient  term  reduced  by  a  factor  of  about  (1  —  a  ),  where  a  =  eft  /E  is  the  ratio  of  well  depth  to 

q'  q  o'  o 

injection  energy. 
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As  grows  beyond  it  begins  to  distort  the  externally-imposed  polyhedral  B  field 
distribution  to  yield  a  configuration  in  which  the  field  is  compressed  outwards  by 
continuing  expansion  of  the  =  1"  surface  (r^).  This  results  in  the  displacement  of 
the  adiabaticity  radius  (r^)  to  larger  radii  at  a  rate  such  that,  with  full  diamagnetic 
currents,  it  will  always  fall  outside  the  ”0  —  1"  radius.  Thus  r^  recedes  outwards  as  r^ 
exceeds  r  . 

X 

Eventually  r^  and  r^  reach  a  critical  radius  rj^  at  which  further  radial  growth  of  r^^ 
becomes  unstable  with  respect  to  r.  This  is  because  the  magnetic  pressure  balance  criterion 
is  satisfied  for  all  radii  beyond  r^^,  once  it  is  reached  at  r^^.  This  radial  position  is  that  at 
which  the  electron  (and  ion)  density  begins  to  increase  rapidly  with  r  >  rj^.  Up  to  r  <  rj^ 

A 

both  ion  and  electron  density  fall  approximately  as  1/r  from  the  core  region  at  radius  r^ 
with  density  n^.  Beyond  r^  both  ion  and  electron  density  increase  rapidly,  and  the  electron 
energy  increases  while  the  ion  energy  decreases  as  r  -» R.  Beyond  the  electron  density 
diverges  markedly  from  the  ion  density  in  order  to  satisfy  Poisson’s  equation  and  produce 
the  desired  negative  potential  wdl.  In  a  general  analysis  of  partide  density  distributions, 
Krall^i  showed  that  the  position  of  <[^>  depends  on  the  the  potential  well  shape.  For 
typical  well  shapes  of  interest,  »  0.83. 

Once  r^^  =  rj^  on  this  model,  the  behavior  of  the  system  "jumps”  from  rj^  to  R,  where 
all  further  increase  in  system  density  buildup  (by  increasing  electron  and  ion  injection 
cunents)  is  stopped  at  the  levflf  relating  drive  current  (I )  and  core  density  (n  )  attained 
at  the  Tj,  =  R  condition.  This  is  (obviously)  because  rj^  can  not  physically  become  greater 
than  R,  even  though  mathematically  this  is  possible,  because  the  B  field  falls  past  r  >  R 
and  the  electron  loss  area  thus  increases  with  inaeasing  r. 


In  short,  the  B  field  can  not  confine  electrons  at  a  density  above  that  at  which  the 
surface  electron  energy  density  exceeds  the  energy  density  of  the  magnetic  field  system. 
Operation  at  this  final  condition  is  pure  WB  mode,  and  is  at  the  point  where  <i^^>  =  1;  all 
mirror  effects  have  vanished  and  the  electron  confinement  is  due  entirely  to  recirculation 
within  the  diamagnetic  sphere  perforated  by  loss  holes  on  the  cusp  axes  whose  radii  are 
proportional  to  the  on-axis  strength  of  the  cusp  B  fields. 
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This  final  <rjj>  =  1  WB  mode  operating  condition  can  be  reached  by  an  alternate 
route.  The  other  extreme  of  startup  and  operation  is  to  avoid  the  MR  mode  entirely,  and 
drive  the  system  from  the  beginning  as  a  high  ^  diamagnetic  device.  This  can  be  done,  in 
principle,  by  starting  with  any  desired  electron  injection  energy  but  with  a  greatly  reduced 
B  field,  so  that  the  /?  ==  1  surface  initially  is  very  far  out  in  the  system,  <8 1-  At  this 
condition  there  will  be  no  mirror  reflection,  because  there  is  no  mirror  field  region;  all  of 
the  field  is  diamagnetically  compressed  to  the  outer  boundary  region,  and  all  electron  losses 
are  governed  by  the  WB  gyro  hole  loss  mechanisms. 


However,  with  small  B  fields  and  large  electron  energy  (as  required  for  the 
production  of  wells  deep  enough  to  yield  ions  with  high  fusion  cross-sections)  the  electron 
recirculation  factor  will  be  very  much  too  small  to  achieve  low  electron  losses  and 
good  power  balance.  In  order  to  reach  high  G.  values  the  confining  B  fidd,  that  sets  the 
size  of  the  electron  gyro  loss  holes  in  the  WB  sphere,  must  be  inaeased  to  a  much  larger 
value.  * 


If  this  is  done  rapidly,  with  no  increase  in  electron  current,  the  /?  =  1  sphere  will  be 
compressed  to  a  radius  smaller  than  r^  s  R,  smd  the  WB  confinement  will  be  reduced  by 
the  increase  in  gyro  hole  loss  area  due  to  the  decreased  gross  B  field  at  smaller  radius,  as 
well  as  by  the  smaller  size  of  the  confining  sphere.  If  this  decrease  is  sufficiently  large,  the 
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system  will  reach  an  unstable  mode  such  that  the  WB  sphere  will  continue  to  shrink  until 
it  stabilizes  with  the  rise  of  MR  mode  confinement  as  r^^  tills  below  rj^.  Now,  the  system 
has  "jumped"  to  operation  in  the  MR  mode,  as  described  above. 

However,  if  the  B  field  is  increased  slowly  (relative  to  electron  lifetime  in  the 
system),  and  in  small  steps,  the  system  will  remain  stable  along  the  =  1  line  and  can 
be  driven  to  higher  and  higher  recirculation  simply  by  increasing  the  B  field,  without 
increasing  the  electron  drive  current.  This  method  avoids  the  limiting  behavior  of  the  MR 
mode,  and  allows  attainment  of  the  final  state  of  maximum  G.  on  the  <r.  >  =  1  line  with 

'  JO  o 

much  less  drive  current  than  is  required  to  make  the  transition  from  MR  to  WB  mode  in 
the  other  startup  scenario. 

In  either  case,  at  this  final  stable  state,  when  =  B.  the  core  and  bulk  densities  of 
electrons  (and  also  of  ions)  are  at  their  maximum  values  for  the  parameters  defining  the 
system  ("e.g.  I^,  n^,  E^,  B^,  R).  Further  increase  in  electron  drive  current  can  not  increase 
system  density.  Since  fusion  power  output  is  proportional  to  the  square  of  the  core  density, 
and  since  this  is  proportional  to  the  surface  density  —  and  hence  to  the  square  of  the  B  field 
—  it  is  clear  that  increasing  B  field  strength  has  a  strong  effect  on  system  fusion  power 
generation.  On  this  simple  argument,  on  the  <z^>  =  1  iine,  the  fusion  power  will  vary  as 
the  fourth  power  of  the  B  field, >3  just  as  for  conventional  magnetic  confinement  machines. 

m.  ANALYTIC  MODELS  OI^LECTRON  MOTION 

In  [I]  the  electron  motion  was  based  on  a  simple  monotonic  potential  well  with 
magnetic  field  and  electron  kinetic  energy  variation  according  to 

B(r)  =  B^<r>“  and  E(r)  =  E^<r>“  (1) 
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where  <i>  s  (r/R).  However,  studies^^*  have  shown  that  the  field  variation 
asymptotically  approaches  the  simple  power  law  form  only  as  the  radius  becomes  smaller 
(approaching  the  origin)  within  the  system.  At  radii  beyond  about  <r>  =  0.5,  the  field 
modulus  exhibits  a  "bumpiness”,  and  its  variation  in  any  plane  section  through  the  center 
of  the  system  and  one  of  the  edge  midpoints  always  shows  a  "rollover"  as  the  edge  of  the 
system  is  approached. 


The  effect  of  this  realistic  departure  from  the  power  law  assumed  in  [I]  is  beneficial 
to  electron  confinement  because  it  gives  larger  B  fields  at  deeper  radii  than  does  the  simple 
power  law  (for  the  same  face  cusp  central  field  strength).  The  simple  model  is  thus  likely 

A* 

to  yield  results  that  are  conservative  (or  pessimistic),  because  of  its  underestimate  of  the 
strength  and  effects  of  the  real  magnetic  fields  in  the  system.  This  can  be  rectified  by 
employing  an  improved  description  of  the  B  field,  which  better  mocks  up  the  actual 
variation  within  the  polyhedron.  Limiting  this,  as  before,  to  radial-only  dependence  (1-D 
can  be  analyzed  in  closed  form;  2-D  can  not)  it  is  found  that  the  edge  effects  can  be 
accounted  for  in  an  approximate  way  by  use  of 


B(r)  =  B^<r>"“f^(r)  =  B^<r>“ 


2 

.(i+<t>"+*). 


(2) 


It  is  obvious  from  this  that  the  field  strength  deep  within  the  system  will  be  twice 
that  previously  used,  for  the  same  value  of  cusp  central-axis  maximum  fidd  strength,  B^. 
As  a  result,  the  electron  gyro  iniii  at  these  inner  radial  positions  will  be  roughly  one-half 
of  those  previously  estimated,  with  concomitant  improvement  in  electron  confinement. 


This  simple  formula  has  the  advantage  that  it  will  still  yidd  analytic  solutions  for 
most  (but  not  all)  of  the  parameters  which  characterize  the  system  and  are  used  in  the 
solution  methods  employed  in  the  modd  of  [I]  (e.g.  the  specific  dimensionless  radii  <r^>i 
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<r.  >,  <r  >,  and  the  tepaxate  mode  xedrculation  latio  terms  G.  and  G.  . 

o  x  jmr  jwD' 

Unfortunately,  the  transit  time  integrals  can  not  be  obtained  in  useful  closed  forms  with 

this  more  realistic  potential  variation.  However,  since  the  transit  time  segments  used  in 

the  calculation  of  the  electron  recirculation  are  divided  by  the  total  transit  time  in  the 

expression  for  the  overall  system  G.  ratio,  the  effect  of  integral  departures  from  exact 

j® 

behavior  are  minimal  for  this  factor,  thus  these  times  can  be  taken  from  the  forms 
described  above  (and  shown  in  Figure  1)  without  serious  error. 

In  the  MR  startup  case,  when  electron  injection  is  started,  the  only  confining 
mechanism  is  mirror— reflection  (MR)  from  the  radius  r^  to  R,  over  a  number  of  transits 
G^.  As  electron  injection  current  is  increased  the  WB  mode  develops  when  r^  >  r^,  and 
the  confinement  parameter  Gj^j^  grows  progressively  larger  than  unity.  Particles  between 
r^  and  r^  simply  transit  through  this  space.  Since  all  particles  that  emerge  from  r^^ 
(escaping  the  WB  region)  enter  the  MR  region  at  r  <  r^,  the  total  number  of  passes  that 
an  electron  will  make  before  loss  is  the  product  of  those  in  each  region. 

To  determine  the  overall  lifetime  and  mean  number  of  passes  made  by  an  electron  in 
the  system  it  is  necessary  to  find  these  parameters  for  each  mode  of  confinement,  and 
weight  them  by  the  fractional  number  of  electrons  participating  in  each  mode.  This 
introduces  a  weighting  factor  dependent  on  <r^>.  With  this  and  the  transit  time  formulae 
summarized  in  Figure  (1),  it  is  possible  to  write*'  the  recirculation  ratios  for  the  two  cases 
when  the  system  has  not  yet  bdfln  to  exhibit  WB  mode  operation  (i.e.  when  r^^  <  r^)  and 
when  the  system  does  show  this  mode  (i.e.  r^  >  r^).  These  formulae  are  given  in 
Figure  (2). 

Using  these  formulae  it  is  possible  to  calculate  G^^  over  a  range  of  the 
system-defining  parameters,  (Z)  and  (W).  Note  that  W  is  defined  by  system  external 


design  criteria,  while  Z  is  determined  by  internal  conditions  reached  in  the  dense  core 
during  system  operation.  It  is  instructive  to  examine  G.  to  see  how  the  MR  and  WB 
contributions  vary  with  these  parameters. 

IV.  ANALYSIS:  THE  GJ  CODE 

The  MR  and  WB  contributions  can  be  calculated  separately  from  their  terms  in  the 
overall  formulae  for  Since  all  of  the  relevant  parametric  terms  in  these  equations  are 
given  only  in  implicit  form,  it  was  necessary  to  construct  a  simple  computer  code  to 
calculate  G.  ,  G.  .  and  G.  .  This  code,  called  the  GJ  code,  was  used  to  determine 
parametric  scaling  of  G^  with  Z  and  W  for  a  wide  range  of  operating  and  design  parameters 
for  the  Polywell^"*  systems  of  interest.  Results  of  these  calculations  are  shown  in  a  series  of 
figures  attached  herewith  (Figures  3-6).  All  calculations  were  made  for  a  truncated  cube 
system  with  an  effective  cusp  number  of  N  s  8,  a  B  field  spatial  index  of  m  s  3,  a  critical 
density  radius  of  <r^>  s  0.83,  loss  radius  factors  of  =  1,  2,  3,  and  for  a  potential  weU 
with  virtual  anode  fractional  height  of  n  s  0.272. 

All  the  figures  show  that  G.  drops  vrith  increasing  Z  (increasing  core  density) 

JQlf 

beyond  the  point  where  r^^  >  r^.  This  is  a  result  of  the  decrease  in  both  MR  mode  time  t^ 
and  in  the  weighting  term  (1— <rj^>),  even  before  <rj^>  reaches  As  <rj^>  moves 

beyond  tbe  drop  is  even  faster,  as  the  diamagnetic  effects  of  induced 

counter-currents  in  the  electr<flftlow  become  larger  and  "push"  <r^>  to  larger  radii. 

Figures  (3a-f)  and  (4)  show  2-D  plots  of  the  total  G.  for  m  =  3  power-law  wells 
with  a  virtual  anode,  over  a  range  of  Z  and  W,  and  for  convergence  ratio  =  lE-2  and 
anode  height  rj  s  0.272.  Figures  (3a-f)  show  a  wide  range  of  the  well  depth  parameter  a^, 
from  *ero  (no  well)  to  0.995.  Figure  (4)  is  a  plot  of  just  the  <rj^>  =  1  line  for  several 
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values  of  the  gyro  loss  radius  parameter  As  expected,  increasing  shifts  the  <rj^>  = 

O 

1  line  to  higher  Z  by  a  factor  of  k^. 

Examination  of  Figures  (3a-f)  shows  that  the  well  depth  affects  the  MR  mode 
drastically  as  becomes  larger  than  about  0.8-0.9,  and  that  the  character  of  the 
parameter  phase  space  changes  for  larger  values  than  0.9.  In  particular,  it  is  found  that  the 
MR  mode  regime  bifurcates  from  the  WB  mode  regime,  and  develops  a  characteristic 
behavior  that  is  antithetical  to  WB  operation.  A  "valley”  appears  between  these  two 
regimes  of  operation  that  subsequent  code  calculations  have  shown  is  very  difficult  to  cross. 

Another  series  of  calculations  was  made  to  test  the  effect  of  anode  height  on  G^. 
These  showed  that  the  variation  of  with  Z  and  W  is  virtually  independent  of  ri  for  all 
T]  >  0.01.  At  very  small  j;  -•  0  the  forms  used  for  the  transit  time  integrals  begin  to  break 
down  and  give  excessive  transit  times;  for  all  realistic  wells  the  dependence  of  G^^  on  anode 
height  is  negligible  for  0  <  t;  <  0.3. 

It  is  useful  to  distinguish  three  regions  of  differing  character  on  these  figures.  In  the 
left  hand  region  below  Z  a  lE13/cm  the  electron  behavior  is  completely  dominated  by 
mirror-reflection  effects.  In  this  area  the  device  is  simply  a  multicusp  mirror  machine, 
operating  with  the  usual  MR  features.  In  the  right  hand  region  above  Z  a  lE16/cm, 
electron  behavior  is  completely  dominated  by  collective  mode  (diamagnetic)  wiffle-ball 
effects,  and  mirror  phenomen#A  this  region  have  essentially  vanished.  Here  and  above  (to 
Z  >  lE18/cm)  lies  the  fusion  reactor  regime.  This  region  is  new  and  unfamiliar  to  the  field 
of  cusp  plasma  research;  it  forms  the  basis  for  the  original  Polywell*”*  concept  and  the  early 
studies^*  of  its  characteristics. 
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In  the  middle  region  where  8E13  <  Z  <  8£15/cm  lies  all  of  the  physics  of  transition 
from  MR  to  WB  operation.  This  region  is  of  experimental  and  theoretical  interest,  because 
of  the  new  and  novel  physics  features  found  here.  However,  this  transition  is  significant  for 
operation  of  these  devices  in  fusion  reactor  regimes  only  if  the  system  is  started  in  the  MR 
mode. 

A  variety  of  calculations  have  been  made  that  show  the  parametric  performance 
across  the  range  from  pure  WB  mode  to  pure  (initial)  MR  mode.  These  are  summarized 
for  one  case  in  Figure  (5)  that  shows  a  three-dimensional  picture  of  the  two-dimensional 
portrayal  of  the  earlier  figures.  From  this  it  is  evident  that  the  device  can  operate  either  as 
a  mirror  machine  OR  as  a  wiffie  ball  machine,  but  that  it  can  NOT  operate  as  both. 

Examination  of  the  calculated  data  shows  that  the  MR  mode  (the  right  hand  area  of 
the  figures)  leads  to  increased  but  only  at  modest  Z,  not  sufhdent  for  power  reactor 
operation.  However,  starting  the  system  as  a  WB  machine  causes  it  to  ride  up  the 
<rj^>  =  1  "ridge"  line  at  the  left  hand  side  of  the  figures,  always  with  an  increasing  Gj. 
Detailed  EKXL  code  calculations  have  been  made,  to  test  the  validity  of  this  GJ  code 
global  mapping,  that  show  the  drive  current  required  to  move  up  either  of  these  startup 
paths. 


The  results  show  that  following  the  MR  mode  at  small  W  leads  to  early  increase  in 
Gy  until  Z  values  around  thos^f  the  -  1  line  are  reached.  Beyond  this  point,  to  force 
the  machine  through  the  transition  from  MR  to  WB  modes,  in  order  to  get  on  the  0—1 
line,  itself,  requires  ever-increasing  electron  drive  currents.  In  contrast,  if  the  device  is 
started  with  a  large  W  (e.g.  by  operation  at  low  B  field),  it  can  be  made  to  reach  the  =  1 
line  at  small  Gj  in  the  WB  mode,  and  then  to  ride  up  the  <r^>  =  1  line  by  slow  reduction 
of  W  (by  slow  increase  of  the  B  field)  to  the  reactor  regime  (ca.  Z  s  lE19/cm)  without 
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significant  increase  in  the  current  required  to  reach  the  =  1  line  in  the  first  place.  This  is 
seen  clearly  in  Figure  7,  where  the  inaeasing  drive  currents  necessary  to  transition  &om 
MR  to  WB  naode  are  plotted  for  a  constant  W.  The  final  current  here  is  18  kA,  where  the 
current  necessary  to  ride  up  the  <r^>  =  1  Hue  is  less  than  1800  A. 

Note  that  the  electron  recirculation  ratio  increases  more  slowly  with  Z,  in  the 
MR  mode,  once  the  Z  value  has  passed  that  for  ^  =  1  (in  the  2— D  projections),  and  that 
the  slope  of  the  G.  vs.  Z  curve,  along  a  line  of  constant  W  first  becomes  negative,  then 
slowly  becomes  positive,  until  it  reaches  a  slope  of  unity.  Below  this  point  the  system  is 
grossly  stable  with  respect  to  increasing  drive  current  input.  However,  beyond  this  point 
the  slope  exceeds  unity,  and  the  system  becomes  grossly  unstable.  From  this  point  up  to 
the  <rjj>  =  1  line  (at  which  point  the  device  can  not  be  driven  further)  the  current  must 
be  reduced,  for  operation  on  the  <r^>  =  1  line  requires  LESS  current  than  to  reach  the 
stable  regions  with  slope  <  1  below  the  unstable  region  under  the  <rj^>  s  1  line.  This  can 
all  be  avoided  by  operation  along  the  /?  =  1  line,  as  described  above. 

Since  the  value  of  G^^  for  MR  operation  drops  with  increasing  Z,  &om  its  initial 
low-Z  value,  as  operation  leaves  the  MR  mode  before  WB  mode  confinement  has  taken 
any  significant  effect  -  it  appears  as  a  "gate”  through  which  machine  operation  must  be 
driven  at  startup.  This  is  true  only  if  started  in  the  MR  mode  and  is  NOT  true  if  started 
in  the  WB  mode.  If  the  device  is  drive-current-4imited,  and  if  the  current  is  insufficient 
(for  the  design  values  of  minimiflF  W  for  the  machine)  to  yield  Z  values  (e.g.  core 
densities)  beyond  this  gate  region,  it  will  not  operate  above  the  much  smaller  Z  value  to 
which  the  limited  current  can  drive  it  in  the  MR  mode. 

The  WB  mode  does  not  begin  until  r^^  >  r^  and  then  rises  rapidly  with  increasing  Z, 
t^j^  and  the  weighting  term  <rj^>.  This  is  quite  clearly  seen  where  the  G^^j^  curves  fall 
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rapidly  to  unity  as  Z  approaches  a  value  of  about  l£14/cm  from  above.  The  cutoff  Z  value 
is  determined  from  the  equation  for  by  setting  this  equal  to  unity.  Thus,  for  all 

A  A 

values  of  Z  <  N(k^)  S  /2r^  the  system  will  be  operating  solely  in  the  MR  mode.  This  has 
the  inverse  consequence  that  WB  mode  physics  can  be  tested  only  if  the  experimental 
system  is  capable  of  being  driven  to  a  core  density  above  this  cutoff  limit; 
n^(r^)^  >  N(kj^)^S^/16itr^.  Typically,  if  r^  =  1  cm  and  N  =  8,  then  n^  >  2ell  1/cm^. 

It  is  very  important  to  note  the  effects  and  characteristics  anticipated  for  this 
device.  In  particular,  a  limited-current  experiment  operating  strictly  in  the  MR  mode  that 
runs  into  the  "gate"  effect  will  demonstrate  poorer  dectron  (and  ion)  confinement  with 
increasing  current  and  core  density.  This  is  to  be  expected  in  this  regime.  Results  of  such 
experiments  can  map  out  the  MR  portion  of  the  3-D  parameter  space  that  describes 
complete  system  operation,  but  can  not  illuminate  the  physics  of  operation  in  the  WB 
regime.  This  regime  is  along  the  <r^>  =  1  ridge  line  in  the  parameter  space  of  the  system, 
as  shown  in  Figure  (5). 

This  final  form  of  the  formalism  outlined  above  (GJ  vers.  1.1),  has  been  built 
into  the  EKXL  code  (vers.  4.1),  so  that  this  code  now  incorporates  fully  correct  physics 
modds  that  describe  system  behavior  across  all  the  possible  range  of  parameter  space. 

Calculations  using  this  version  4.1  of  the  EKXL  code  have  been  made  paramet¬ 
rically  for  a  variety  of  cases  at H^h  density  ("reactor")  conditions,  and  for  a  wide  range  of 
modds  for  SCIF  experimental  systems  and  operations.  Results  of  the  SCIF  computations 
are  discussed  in  a  separate  EMC2  technical  note^^*.  Results  of  the  reactor  case  studies  are 
given  in  a  forthcoming  EMC2  technical  note.  In  all  of  these  the  dfective  "loss  radius"  has 
been  taken  to  be  in  the  range  2  <  <  3,  as  derived  from  dementaiy  considerations^^  of 
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diamagnetic  behavior  in  the  system.  Knowledge  derived  from  this  work  has  been  used  to 
analyze  possibilities  for  startup  of  SCIF  with  neutral  background  gas;  this  is  reported  in 
another  technical  note.33 


LIST  OF  FIGURES 


figure  1. 

Summary  of  formulae  used  for  total  and  partial  (regmented)  electron  transit  time  calculation  in 
virtual  anode  wells  with  "rollover"  B  field. 


Figure  2. 


Summary  of  equations  for  total  electron  recirculation  ratio,  and  under  mirror  reflection  (MR)  and 
wiffle— ball  (WB)  operating  modes,  and  formulae  used  for  critical  system  radii  (gyro,  adiabaticity, 
transition  and  magnetic  pressure  balance  radii)  in  "rollover"  wells. 


Figure  3.  a— f 

Effect  of  well  depth,  a  ,  on  G.  vs.  Z  plots  for  various  values  of  W.  The  <r,  >  si  1  line  is 

q  J  b 

indicated  in  each  chart  by  the  dashed  line. 


Figure  4. 


Effect  of  gyro  "hole"  radius,  k^,  on  Gj  vs.  Z  plots  for  various  values  of  W.  a)  The  =  1 


line  for  =  1,  2,  3.  The  magnitude  of  W  at  Gj  s  1  is  indicated. 


Figure  5. 


3— D  plot  of  2— D  projections  with  k^— 2. 


Figure  6.  &,b 

Increasing  drive  current  necessary  to  transition  from  MR  to  WB  mode  in  a  system  with  a  fixed 
W.  Currents  given  are  in  kA. 
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Summary  of  equations  for  total  dectron  recirculation  ratio,  and  under  mirror 
redectioiv  (MIQ  and  wiffle-ball  (WB)  operating  modes,  and  formulae  used  for 
critical  system  radii  (^ro,  adiabaticity,  transition  and  magnetic  pressure  balance 
radii)  in ''rollover”  wdls. 
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Effect  of  weU  depth,  on  G.  vs.  Z  plots  for  various  values  of  W.  The  <r, 
line  is  indicated  in  each  chart  by  the  dashed  line. 


3_D  plot  of  2-D  projections 
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Figure  6  a. 
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Figure  6  b. 
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Figure  6.  a,b 

Increasing  drive  current  necessary  to  transition  from  MR  to  WB  mode  in  a  system 
with  a  fix«i  W.  Cunents  given  are  in  kA. 


